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ABSTRACT 

We present a compilation of properties of the 105,783 quasars in the SDSS Data Release 7 (DR7) quasar 
catalog. In this product, we compile continuum and emission line measurements around the Ha, H/?, Mgii 
and Civ regions, as well as other quantities such as radio properties, and flags indicating broad absorption line 
quasars (BALQSOs), disk emitters, etc. We also compile virial black hole mass estimates based on various 
calibrations. For the fiducial virial mass estimates we use the Vestergaard & Peterson (VP06) calibrations for 
H/3 and Civ, and our own calibration for Mgii which matches the VP06 H/? masses on average. We describe 
the construction of this catalog, and discuss its limitations. The catalog and its future updates will be made 
publicly available online. 

Subject headings: black hole physics — galaxies: active — quasars: general — surveys 



1. INTRODUCTION 

In recent years, studies of quasars and active galactic 
nuclei (AGNs) have been greatly facilitated by dedicated 
large-scale wide and deep field surveys in different bands, 
most notably by optica l surveys such as the Sloan Digi- 
tal S ky Survey (SDSS |^AeLaLl|2QQP^ and the 2QZ sur- 
vey (I Croorii et alj|2004|) . Indeed, the growing body of data 
has revolutionized the study of quasars and AGNs. Large, 
homogeneous data sets allow detailed investigations of the 
phenomenological properties of quasars and AGNs, offer- 
ing new insights into the central engine powering these ob- 
jects and their connections to their host galaxies, especially 
when combined with multi- wavelength coverage. At the 
same time, it has become important to fit the quasar/ AGN 
population into its cosmological context, i.e., how the su- 
permassive black hole (SMBH) population evolves across 
cosmic time. These data have led to a coherent picture 
of the cosmic evolution of the S MBH population within 
the c o ncordance ACDM para d igm ("Kauffmann & Haehnelti 
2(M IWyiflie & Loebl llopl iHopkins et al. 2006, 2008; 



Shankar et alJl2009l:IShenll2009h . where the key observational 



components include quasar clustering, the luminosity func- 
tion (LF), the BH mass function, and the correlations between 
BHs and their host properties. Increasingly larger data sets are 
offering unique opportunities to measure these properties with 
unprecedented precision. 

In an earlier study of the virial BH mass and Eddington 
ratio distributions of quasars, we measured spectral prop- 
erties for the SPSS Data Release 5 (DR 5) quasar catalog 
jSchneider et alj 120071: IShen et alj l20()8bl) . We hereby ex- 
tend t his exercise to the DR7 quasar catalog jSchneider et alj 
I20T0I) . We now include a more complete compilation than be- 
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fore of quantities from our spectral fits. Our measurements are 
more sophisticated than the SDSS pipeline outputs in many 
ways, and are hence of practical value. We describe the par- 
ent quasar sample in ^ and the spectral measurements and 
the catalog format in |3] We discuss possible applications of 
our measurements in ^ Throughout this paper we use cos- 
mological parameters J^a = 0.7, flo = 0.3 and h = 0.7. 

2. THE sample 

The SDSS use s a dedicated 2.5-m wide-field telescope 
jGunn et al.! '2006') with a drift-scan camera with 30 2048 x 
2048 CCDs (Gunnetal. 1998) to imag^ the sky in five 
broad bands {ugriz; iFukugita et alj 119961) . The imag- 
ing dat a are taken on d ark photometric nights of good 
seeing jHogg et alJ l2001h. are cal i brated photometr ically 
("Smith et al."2002';'Ivezic et al ."2004; 'Tucker et al."2006) and 
astrometrically (Pieretal. 2003), and object parameters are 
measured (Lupton et al. 2001; Stou ghton et al.ll2002l) . Quasar 
candidates ( Richards et al.. 2002a) for follow-up spectroscopy 
are selected from the imaging dat a using their colors, a nd are 
arranged in spectroscopic plates (iBlanton et al.l 120031) to be 
observed with a pair of fiber-fed double spectrographs. 

Our parent sample is the latest co mpilation of the spectro- 
scopic quasar catalog ( Schneider et al. 2010) from SDSS DR7 
(lAbazajian et al J 120091) . This sample contains 105,783 bona 
fide quasars brighter than M,- = -22.0 and have at least one 
broad emission line with FWHM larger than 1000 km s"' or 
have interesting/complex absorption features. About half of 
these objects are selected uniformly using the final quasar tar - 
get selection algorithm described in [Richards et alj (l2002al) . 
with the remaining objects selected via early versions of 
the target se lection algorithm or various serendipitous algo- 
rithms (see Schneider et al. 2010), whose selection complete- 
ness cannot be readily quantified. For statistical studies such 
as quasar clustering and the LF, one should use the uniformly 
selected quasar sample. Fig. [T] shows the distribution of the 
105,783 quasars in the redshift-luminosity plane. 

The reduced one-dimensional (ID) spectral data used in 
this study are available through the SDSS Data Archive 
Server^ (DAS). The spectral resolution is - 1850-2200, 

^ http://das.sdss.org/spectro/ 
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and the ID spectra are stored in vacuum wavelength, with 
a pixel scale of 10"^ in log-wavelength, which corresponds 
to ^ 69 kms"'. Since the 6t h SDSS data release (DR6, 
lAdelman-McCarthv et al.ll2008h . the spectral flux calibration 
is scaled to the point spread function (PSF) magnitudes of 
standard stars, therefore there is no longer need for a fiber- 
to-PSF conversion for the spectral flux (Shen et al. 2008 also 
used the PSF spectral flux calibration). Throughout the paper, 
we refer to the signal-to-noise ratio per pixel as S/N. 

To include radio prope rties, we match th e DR7 quasar 
catalog with the FIRST jWhite etalj|1997h catalog*^ with 
a matching radius of 3 0" and estimate th e radio loudness 
R = /6cm//2500 following Ijiang et aP (l2007h . where f(,cm and 
/2500 are the flux density (/j^) at rest-frame 6 cm and 2500 
A, respectively. For the quasars that have only one FIRST 
source within 30" we match them again to the FIRST catalog 
with a matching radius 5" and classify the matched ones as 
core-dominant radio quasars. Those quasars that have multi- 
ple FIRST source matches within 30" are classified as lobe- 
dominated. The rest-frame 6 cm flux density is determined 
from the FIRST integrated flux density at 20 cm assuming 
a power-law slope of a^, = -0.5; the rest-frame 2500 A flux 
density is determined from the power-law continuum fit to 
the spectrum as described in ^ For lobe-dominated radio 
quasars we use all the matched FIRST sources to compute the 
radio flux density. We note that we may have missed some 
double-lobed radio quasars with lo be diameter larger th an 1 ' . 

To flag BALQSOs, we use the i Gibson et al.1 (l2009h DR5 
BALQSO catalog to set the Civ and Mgii BALQSO flags (us- 
ing their "BIO" flags). We also visually inspected all the post- 
DR5 quasars with redshift z > 1 .45 to identify obvious Civ 
BALQSOs (we may have missed some weak BALQSOs). We 
did not perform a systematic search for low-ionization Mgii 
BALQSOs because of the large number of quasars with Mgii 
coverage and the much rarer occurrence of Mgii BALQSOs. 
Although we report serendipitously identified Mgii BALQ- 
SOs, the completeness of these objects is low. We identified a 
total of 6214 BALQSOs in the DR7 quasar catalog. 

There are also subclasses of quasars which show interest- 
ing spectral features in their broad or narrow emission line 
profiles. Some quasars show a double-peaked or asymmetric 
broad Balmer line profile, which is commonly interpreted as 
arising from a relativ istic accretio n disk arou nd the black hole 
(disk emitters, e.g.. IChen et al.l [1989; Erac leous & HalpernI 
[1994; Strateva et al. 2003) although alternative interpretations 
exist for some of these objects, such as a binary SMBH (e.g., 
lGaskellll2010l and references therein). Some quasars show 
dou ble-peak ed narro w lines (such as [O iii] AA4959,5007 , 
e.g.,[Liuetal. 2010al [SmMTetH] [20091 : IWang et alJ[2009ah . 
which could be due to either narrow line region kinematics or 
a merging AGN pair (e.g., Liu et al. 2010b; Shen_et ak 201 1). 
We have visually inspected all of the z < 0.89 quasars in the 
DR7 catalog (i.e., those with H/3 coverage) and flagged such 
objects. 

3. SPECTRAL MEASUREMENTS 

We are primarily interested in the broad Ha, H/3, Mgii, 
and Civ emission lines because these are the most fre- 
quently studied lines that are available over a wide range 
of redshifts. More importantly, they have been cali- 
brated as virial black hole (BH) mass estimators (e.g., 

* The version of the FIRST source catalog and the coverage maps used are 
as of July 16, 2008 l http://sundog. stsci.edu/first/catalogs/readme. html I. 



-30 




12 3 4 5 

Redshift 



Fig. I. — The distribution of DR7 quasars in luminosity-redshift space, 
where the left and bottom panels show the luminosity and redshift histograms. 
Luminosity is indicated using the (continuum and emission line) ^T-coirected, 
(■-band absolute magnitude, Mi(z = 2), normalized at z = 2 (Richards et al. 
2006a). Th e bla ck dots are uniformly- selected quasars (see Richards et al. 
2002a and j|3.8t and the gray dots are quasars selected by a variety of ear- 
lier algorithms or serendipitous selections. The red dots are selected by the 
QSO_HiZ uniform selection (e.g., Richards et al. 2002a) but with ; > 19.1 
and at z < 2.9, and should be removed in constructing homogeneous quasar 
samples. The cyan lines show the con'esponding M,(z = 2) for = 19. 1 
(z < 2.9) and = 20.2 (z > 2.9) respectively, and the gray dashed line shows 
the equivalent for i = 15 (the bright limit for SDSS quasar targets). The non- 
uniformly selected quasars (gray dots) are targeted to fainter luminosities than 
are the uniformly selected quasars. 

Vestergaardl2002HMcLure & Jarvisll2002l;lMcLure & Dunlod 
'2004"; 'Greene & Ho" '2005b'; 'Vestergaard & Peterson' '20061 
McGifl et al. 2008; Vestergaard & Osmer 2009; Wang et aQ 
2009b). 

There are numerous studies of the statistical emission line 
properties of quasars relying either on dir ect measurements, 
or ori spectral fits of the line profile (e.g., Boroson & Greed 
119921; Marziani et al. 1996; McLure & Jarvis 200 J 
Richards et al. 2002b; McLure & Dunlop 2004; Bachev et at| 
2004; ,Dieti-ich & Hamann 2004; Baskin & L aor ,20051 
Kollmeier et'aiT[2 006; Fine et al. 2006, 2008; Bo nning et al. 
'2007; 'Salviander et al. 2007; Sulentic et al. 2007; Shen et aL 
2008a.b; Hu et al. 2008a b; Zamfiretal. 2010; Wang et al" 
l2009bt IWuetal.ll2009t iDong et aT.ll2009bllal) . For the same 
set of data, different studies sometimes report different 
results for certain measured quantities due to the different 
line-measurement techniques used in these studies. Which 
method is preferred, however, depends on the nature of the 
problem under study. A classic example is measuring the 
full-width-at-half-maximum (FWHM) in estimating the BH 
mass using virial estimators, where the usual complications 
are: 1) how to subtract the continuum underneath the line; 2) 
how to treat the narrow line component (especially for Mgii 
and Civ); 3) how to measure the broad line profile, especially 
in the presence of noise and absorption. These choices 
crucially depend on the particular virial estimator calibrations 
used; in particular, one must use a similar FWHM definition 
as was used in the virial estimator calibration, and new meth- 
ods of FWHM measurements must be re-calibrated either 
against reverberation mapping (RM) masses or internally 
between different line estimators. On the other hand, different 
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line-measurement methods have different sensitivities to the 
quality of the spectra (spectral resolution and S/N), which 
introduce systematics when switching fro m high-quality to 
low-quality data (e.g., iDennev et al.ll2009l) . It is beyond our 
scope to fully settle these issues within the current study. 

We remove the effects of Gal actic e xtinction in the SDSS 
spectra using the Schlegel et a l.] (Il998h map a nd a Milky Way 
extinction curve from C ardelUet al.1 (Il989h with Ry = 3.1, 
and shift the spectra to restframe using the cataloged red- 
shift as the systemic redshift^. For each line, we fit a lo- 
cal p ower-law continuum (/> =AA°^) plus an iron tem- 
plate (iBoroson & Greed 119921: IVestergaard & WilkesI 120011: 
ISalviander et al.l l2007 l) to the wavelength range around the 
line that is not contaminated by the broad line emission. Dur- 
ing the continuum+iron fitting we simultaneously fit five pa- 
rameters: the normalization A and slope ax of the power- 
law continuum, and the normalization Ape, line broadening 
(Tpe and velocity offset vpe relative to the systemic redshift 
for the iron template fit. Because of the moderate spec- 
tral quality of SDSS spectra (median S/N< 10) (Tpe and vpe 
are often poorly constrained; nevertheless the iron fit gives 
a reasonably good estimate of the iron flux to be subtracted 
off. The continuum+iron fit is then subtracted from the spec- 
trum, and the resulting line spectrum is modelled by various 
functions. In the case of Ha and H/3 the adjacent narrow 
emission lines, e.g., [Oiii] AA4959,5007, [Nil] AA6548,6584, 
[S ii] AA67 17,6731, are also fit simultaneously. Below we de- 
scribe the detailed fitting procedures for the four broad lines. 

3.1. Ha 

For Ha we use the optical iron template from 
iBoroson & GreenI (|1992|) . and we fit for objects with z < 0.39. 
The continuum+iron fitting windows are [6000,6250] A and 
[6800,7000] A. 

For Ha line fitting, we fit the wavelength range [6400,6800] 
A. The naiTow components of Ha, [Nil] AA6548,6584, 
[Sii] AA6717,6731 are each fit with a single Gaussian. Their 
velocity offsets from the systemic redshift and line widths 
are constrained to be the same, and the relative flux ratio 
of the two [Nil] components is fixed to 2.96. We impose 
an upper limit on the narrow line FWHM < 1200 kms"' 
(e.g.,lHao et al. 2005). The broad Ha component is modelled 
in two different ways: a) a single Gaussian with a FWHM 
> 1200 kms"'; b) multiple Gaussians with up to three Gaus- 
sians, each with a FWHM > 1 200 km s~' . The second method 
yields similar results to the fits with a truncated Gaussian- 
Hermite function (e.g., van der Marel & Franx 1993). During 
the fitting, all lines are restricted to be emission lines (i.e., 
positive flux). 

3.2. H/S 

For H/j we use the optical iron template from 
IBoroson & GreenI (Il992l) . and we fit for objects with z < 0.89. 
The continuum+iron fitting windows are [4435,4700] A and 
[5100,5535] A. For the H/3 line fitting, we follow a similar 
procedure as Ha to fit for H^ and [Oiii] AA4959,5007, where 
the line fitting wavelength range is [4700,5100] A. Since 
the [O III] AA4959, 5007 lines frequ ently show asymmetric 
blue wings (e.g.. iHeckman et al.llT98ll: iGreene & Holl2005al: 
iKomossa etal.ll2008h and sometimes even more dramatic 

' IHewett & Wlia OoTQ") provided improved redshifts for SDSS quasars. 
However, this subtlety has negligible effects on our spectral fits. 



double-peaked profiles (e .g., iLiu et al.l 1201 Oat ISmith et aTl 
I2OO9I: IWang et al.l l2009ah . we model each of the narrow 
[O III] AA4959,5007 fines with two Gaussians, one for the 
core and the other for the blue wing. The flux ratio of the 
[O III] doublet is not fixed during the fit, but we found that the 
fitting results show good agreement with the theoretical ratio 
of about 3. The velocity offset and FWHM of the narrow 
H/3 line ai-e tied to those of the core [O iii] AA4959,5007 
components'", and we impose an upper limit of 1200 kms"' 
on the narrow line FWHM. As in the Ha case, the broad H/3 
component is modelled in two ways: either by a single Gaus- 
sian, or by multiple Gaussians with up to three Gaussians, 
each with a FWHM > 1200 kms"' . 

The single Gaussian fit to the br oad component is e ssen- 
tially the same to the procedure in Shen et al. (2008b), and 
is somewhat similar to the procedure in McLure & DunloJ 
(2004)". However, in many objects the broad Ha/H/3 com- 
ponent cannot be fit perfectly with a single Gaussian; and 
FWHMs from the single Gaussian fits are systematically 
large r by ~ 0. 1 dex than those from the multiple Gaussian fits 
(e.g.. IShen et al.ll2008bh . The additional multiple Gaussian 
fits for the broad Ha/H/3 component provide a better fit to the 
overall broad line profile, and the FWHM measured from the 
model flux can be used in customized virial calibrations. It 
is unclear, however, which FWHM is a better surrogate for 
the true virial velocity, that is, the one that yields the smallest 
scatter in the calibration against RM black hole masses. 

3.3. Mgu 

For Mgii we u s e the UV iron template from 
IVestergaard & WilkesI (12001 1) , and we fit for objects with 
0.35 < z < 2.25. The continuum+iron fitting windows are 
[2200,2700] A and [2900,3090] A. We then subtract the 
pseudo-continuum from the spectrum, and fit for the Mgii 
line over the [2700,2900] A wavelength range, with a single 
Gaussian (with FWHM < 1200 kms"') for the narrow Mgii 
component, and for the broad Mgii component with: 1) 
a single Gaussian; 2) multiple Gaussians with up to three 
Gaussians. Again, the multiple-Gaussian fits often provide 
a better fit to the overall broad Mgii profile; but we have 
retained the FWHMs from a single Gaussian fit in orde r to use 
the Mgii virial m ass calibrations in iMcLure & JarvisI (l2002l) 
and iMcLure & D unlop (2004). S ome Mgii virial estimator 
calibrations (e.g., McLure & Jarvis 20021; iMcLure & DunlopI 
,2004; Wang et al. 2009b) d o subtract a narro w Mgii com- 
ponent while others (e.g., IVestergaard & Os mer 20 0^ do 
not. T o utilize the Mgii calibration in IVestergaard & Osmerl 
(l2009l) . we also measure the FWHMs from the broad+narrow 
Mgii fits (with multiple Gaussians for the broad component), 
where any Gaussian component having flux less than 5% of 
the total line flux is rejected when computing the FWHM — 
this step is to eliminate artificial noise spikes which can bias 
the FWHM measurements. During our fitting, we mask out 

In some objects the narrow H/3 component might have different width 
and velocity offset from those of the core [O III] AA4959,5007 components. 
Usually in such cases our procedure still provides a reasonable approximation 
to subtract the narrow H/3 component; however, under rare circumstances 
where the naiTow H/3 component is not obvious this may lead to a biased 
narrow line subtraction. Nevertheless, fitting the narrow H/3 line without this 
constraint would significantly degrade the reliability of narrow line subtrac- 
tion. 

' ' In addition to Gaussian profiles. [McLure & Dunloij j2004l) also tried to 
fit the broad/narrow component with a single Lorentzian, but this does not 
change the measured broad FWHM significantly. 
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Fig. 2. — Two examples of Mgll lines which show narrow line components. 
The spectra are plotted as black lines and the gray lines show the errors per 
pixel. The dashed vertical lines mark the locations of the Mgll AA2796,2803 
doublet. 

3(7 outliers below the 20-pixel boxcar-smoothed spectrum to 
reduce the effects of narrow absorption troughs. 

Unlike the cases of Ha and H/3, it is somewhat ambigu- 
ous whether it is necessary to subtract a narrow line com- 
ponent for Mgii and if so, how to do it. On one hand, 
for some objects, such as SDSSJ002209.95H-001629.3 and 
SDSSJ130101.94H-130227.1 (e.g., Fig.Ell, the spectral quality 
is sufficient to see the bifurcation of the Mgii doublet around 
the peak. The locations of the two peaks indicate that they are 
associated with the Mgii AA2796,2803 doublet, and the fact 
that they are resolved means that the FWHM of each com- 
ponent is < 750 kms"', hence they are most likely associ- 
ated with the narrow line region. On the other hand, such 
cases are rare and most SDSS spectra do not have adequate 
S/N to unambiguously locate the narrow Mgii doublet. As- 
sociated narrow Mgii absorption troughs can further compli- 
cate the situation by mimicing two peaks. Hence although 
our approach of fitting a single Gaussian to the narrow Mgii 
component is not perfect, it nevertheless accounts for some 
narrow Mgii contamination. Fig.[3]compares our broad Mgii 
FWHM measurements with those from Wang et al. (2009y) 
for the objects in both studies. Although we have used a dif- 
ferent approach, our results are consistent with theirs, with a 
mean offset ^ 0.05 dex. This systematic offset between our 
results and theirs is caused by the fact that they are treating 
the broad Mgii line as a doublet as well, while we (and most 
studies) are treating the broad Mgii line as a single compo- 
nent. 

3.4. Civ 

For Civ we fit for objects with 1.5 < z < 4.95. Iron emis- 
sion is generally weak for Civ, and most of our objects do not 



have the spectral quality sufficient f or a reliable iron template 
subtraction (e.g., IShen et aTll2008bl) . The power-law contin- 
uum fitting windows are: [1445,1465] A and [1700,1705] A. 
We found fitting Civ with iron subtraction does not change 
the fitted Civ FWHM significantly, but does systematically 
reduce the Civ EW by ~ 0.05 dex because the iron flux un- 
der the wings of the Civ line is accounted for. At the same 
time, fitting iron emission increases the uncertainty in the fit- 
ted continuum slope and normalization, due to imperfect sub- 
traction of the iron flux. Therefore we report our Civ mea- 
surements without the iron template fits, and emphasize that 
the Civ EWs may be overestimated by ^ 0.05 dex on average. 

The continuum subtracted line emission within [1 500,1600] 
A was fitted with three Gaussians (e.g., IShen et al. 2008b), 
and we measure the line FWHM from the model fit. To re- 
duce the effects of noise spikes, we reject any Gaussian com- 
ponent having flux less than 5% of the total model flux when 
computing the FWHM. However, unlike some attempts in 
the literature (e.g., |B achev et al. 2004; Baskin & Laor 2005] 
ISulentic et a f 2007; ZamfiretaLl|2010), we do not subtract a 
narrow Civ component because: 1) it is still debatable if a 
strong narrow Civ component exists for most quasars, or if it 
is feasible to do such a subtraction; 2) existing Civ virial es- 
timator s are calibrated using the FWH Ms from the entire Civ 
profile (IVestergaard & Petersonll2006 ). 

Many Civ lines are affected by narrow or broad absorp- 
tion features. To reduce the effects of such absorption on the 
Civ fits, we mask out 3<t outliers below the 20-pixel boxcar- 
smoothed spectrum during our fits (to remedy for narrow ab- 
sorption features); we also perform a second fit excluding pix- 
els below 3(7 of the first model fit, and replace the first one 
if statistically justified (to account for broad absorption fea- 
tures). We found these recipes can alleviate the impact of 
narrow or moderate absorption features, but the improvement 
is marginal for objects severely affected by broad absorption. 

3.5. Reliability of spectral fits and error estimation 

Our spectral fits were performed in an automatic fashion. 
Upon visual inspection of the fitting results we are confident 
that the vast majority (> 95%) of the fits to high S/N spectra 
were successful, and comparisons with independent fits by 
others also show good agreement. However, the reliability of 
our spectral fits drops rapidly for low-quality spectra. Fig. 
m shows the distributions of the median S/N per pixel around 
the line-fitting region for objects that have line measurements, 
for Ha, H/3, Mgii and Civ respectively. Although the bulk of 
objects have median S/N> 5 for the line-fitting regions, there 
are many objects that have lower median S/N, especially for 
Civ at high redshift. The effects of S/N on the measurements 
depend on both the properties of the lines (i.e., line profile, 
line strength, degree of absorption features, etc), and the line- 
fitting technique itself (i.e., what functional form was used, 
how to deal with absorption troughs, etc). 

To investigate the impact of S/N on our fitting parameters 
we ran a series of Monte Carlo simulations. We select repre- 
sentative real spectra with high S/N, then degrade the spectra 
by adding Gaussian noise and measure the line properties us- 
ing the same line-fitting routine. For each line (H/3, Mgii, or 
Civ), we study several objects with various line shapes and 
EWs. We simulate 500 trials for each S/N level and take the 
median and the 68% range as the measurement result and its 
error. 

Figs. |5]|8] show several examples of our investigations for 
Ha, H/3, Mgii and Civ respectively. As expected, decreasing 
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Fig. 3. — Comparison of our Mgll FWHMs with tliose from lWang etaj] <20Q9W for the same objects. The left panel shows a scatter plot, where the solid Hne 
is the unity relatio n. Th e right panel shows a histogram of the ratio between the two values. Our broad Mgll FWHM values are systematically larger by ~ 0.05 
dex than those in lWang et alj l2009bi) . mainly caused by the fact that they fit the broad Mgll as a doublet while we did not. 

noise to the original spectrum using the reported flux den- 
sity errors, and fit for those mock spectra with the same fitting 
routines. We estimate the measurement uncertainties from the 
68% range (centered on the median) of the distributions of fit- 
ting resuhs of the 50 trials. This approach gives more reason- 
able error estimation than using the statistical errors resulting 
from the ^ fits, in the sense that it not only takes into ac- 
count the spectral S/N but also the ambiguity of subtracting a 
narrow line component in many cases (for Ha, H/3 and Mgii). 
We estimate the spectral measurement uncertainties for all the 
105,783 quasars in our sample in this way. 

Given the automatic nature and specifics of our fitting 
recipe, there will undoubtedly be bad fits for noisy spectra or 
peculiar objects. We recommend using the reported measure- 
ment errors to remove suspicious measurements, e.g., Err> 
AND E r r < some threshold. Our fitting recipe was optimized 
for the vast majority of quasars in our sample, and it may fail 
badly for rare objects with peculiar continuum and emission 
line properties. These objects include severe BALQSOs, disk 
emitters, and objects such as J094215.12H-090015.8, which 
has extremely broad lines that exceed our line fitting range. 
One should pay attention to our fitting ranges and these spe- 
cial objects upon usage of the cataloged quantities. For var- 
ious reasons, one may want to check the quality assessment 
(QA) plot for individual fits to make sure that the fitting re- 
sults are robust. Such QA plots are provided along with this 
catalog. 

Finally, we note that we imposed a FWHM limit of 
1200 kms"' for the narrow line components when fitting Ha, 
H/3 and Mgii. This choice is motivated by the results in 
iHao et all (|2005|) . but is still somewhat arbitrary; narrow lines 
broader than this threshold are not unusual. On the other hand, 
if a Gaussian component with FWHM< 1200kms"' can be 
fit to Ha, H/3 and Mgii, it will be considered as a narrow 
line component and subtracted off. For these reasons, we 
urge caution regarding objects for which the measured narrow 
line FWHM reaches the 1200 kms"' limit (^ a few percent). 
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Fig. 4. — Distributions of the median S/N per pixel around the line-fitting 
regions for objects with line measurements, for Ha, H/3, Mgll and CIV re- 
spectively, where the underlying area of each distribution is normalized at 
unity. 

the S/N increases measurement scatter. In all cases the fitted 
continuum is unbiased as S/N decreases. The FWHMs and 
EWs are biased by less than ±20% for high-EW objects as 
S/N is reduced to as low as ^ 3. For low-EW objects, the 
FWHMs and EWs ai-e biased low/high by > 20% for S/N< 5. 
Since the median EWs for the three lines are > 30 A (see 
^4.11 ). we expect that the measurements for most objects are 
unbiased to within ±20% down to S/N^ 3. But for many 
purposes, it would be more conservative to impose a cut at 
S/N> 5 for reliable measurements. 

To estimate the uncertainties in the measured quantities in 
our fits, we generate 50 mock spectra by adding Gaussian 
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Fig. 5. — Effects of S/N on the line measurements for Ha for three rep- 
resentative examples. For each object we show the actual spectrum (black 
hne) and the best-fit model (red line) in the upper panel. The lower three 
panels show the ratios of the values measured from the degraded spectra to 
those measured from the original spectrum, as functions of S/N; black dots 
are median values and the error bars indicate the 68% quantile. 

Upon visual inspection of the QA plots for these objects, even 
though the FWHM parameter reaches the limit, most of the 
fits still yield reasonable measurements of total Une flux, line 
centroid and narrow line subtraction. 

3.6. Host galaxy contamination 
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Fig. 6. — Effects of S/N on the line measurements for H/3 for three rep- 
resentative examples. For each object we show the actual spectrum (black 
line) and the best-fit model (red line) in the upper panel. The lower three 
panels show the ratios of the values measured from the degraded spectra to 
those measured from the original spectrum, as functions of S/N; black dots 
are median values and the error bars indicate the 68% quantile. 

For the vast majority of objects in our catalog with z > 0.5, 
host galaxy contamination is negligible. However, for the 
z < 0.5 low-luminosity quasars in our sample, the continuum 
luminosity at restframe 5100A may be contaminated by light 
from the host galaxies. Unfortunately the spectral quality of 
the majority of individual objects does not allow a reliable 
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Fig. 7. — Effects of S/N on the line measurements for Mgll for three 
representative examples. For each object we show the actual spectrum (black 
hne) and the best-fit model (red line) in the upper panel. The lower three 
panels show the ratios of the values measured from the degraded spectra to 
those measured from the original spectrum, as functions of S/N; black dots 
are median values and the error bars indicate the 68% quantile. Note that our 
continuum+iron fit does not account for the [Ne iv]/Fe iii and [O ii] emission 
around 2400 A-2480 A. 
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Fig. 8. — Effects of S/N on the line measurements for CIV for three rep- 
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Fig. 9. — Composite spectra for objects binned in logLjioo, normalized at 
3000A. The flux gradually flattens at long wavelength due to increasing host 
contamination towards fainter quasar luminosities, accompanied by the in- 
creasing prominence of stellar absorption features and narrow line emission. 
The inset shows the fractional host contamination at 5 lOOA assuming that the 
highest luminosity bin (logLsjoo = 45.5) is not affected by host emission and 
that the intrinsic AGN power-law continuum slope does not change over the 
luminosity range considered. The dashed line in the inset is a polynomial fit 
(Eqn.[TJ. 

galaxy continuum subtraction. Here we estimate the effects 
of host contamination with stacked spectra. 

We take all quasars with measurable rest frame 5100A 
continuum luminosity, L5100, and bin them on a grid of 
Al ogLsioo = 0.1 for logCLsion /ergs"^) = 44.1-45.5. Follow- 
ing I Vanden Berk et al] (1200 Ih . we generate geometric mean 
composite spectra for objects in each luminosity bin. The 
composite spectra are shown in Fig. |9] where the flux grad- 
ually flattens at long wavelengths due to increasing host con- 
tamination at fainter luminosities. This trend is accompanied 
by the increasing prominence of stellar absorption features 
and narrow line emission towards fainter luminosities. The 
inset shows the fractional host contamination at 5100A as- 
suming that the highest luminosity bin (logLsloo = 45.5) is 
not affected by the host and that the intrinsic AGN power- 
law continuum slope does not change over the luminosity 
range considered. The host contamination is substantial at 
logLsioo < 44.5, and becomes negligible towards higher lu- 
minosities. The median value of logLsloo for quasars in this 
low-redshift sample is ~ 44.6, and therefore the host contam- 
ination on average is ^ 15%, which leads to a ^ 0.06 dex 
overestimation of the 5 100 A continuum luminosity and thus 
^ 0.03 dex overestimation of the H/3-based virial masses for 
the median object. While we do not correct the measured 
5100A continuum luminosity (and other quantities depending 
on it) in the catalog, we provide an empirical fitting formula of 
the average host contamination based on the stacked spectra 
(dashed line in the inset of Fig. |9l): 



.^SlOO.host 



: 0.8052- 1.5502x+0.9121x^-0.1577x' 



(1) 



Lsm.Qso 

for x + 44 = log(L5ioo,totai/ergs~') < 45.053; no correction is 
needed for luminosities above this value. 

We suspect that host contamination is largely responsible 
for the apparent anti-correlation between L5100 and spectral 



slope a\, and the "negative" Baldwin effect (e.g., iBaldwinI 
Il977h between t he b road H/3 EW and L5100 below L5100 ~ 
10^5 ergs ' (see ^4. 11 1. 



3.7. Virial BH masses 



It has become common practice to estimate quasar/ AGN 
BH masses based on single-epoch spectra (hereafter virial 
mass in short). This approach assumes that the broad line re- 
gion (BLR) is virialized, the continuum luminosity'^ is used 
as a proxy for the BLR radius, and the broad line width 
(FWHM or line dispersion) is used as a proxy for the virial 
velocity. The virial mass estimate can be expressed as: 

Mq J ^VlO ergs-iy ^Vkms-iy' 

(2) 

where the coefficients a and b are empirically calibrated 
against local AGNs with RM masses or internally among 
different lines. H/3, Mgii, Civ, and their corresponding 
continuum luminosities are all frequently adopted in such 
virial calibrations. Although it is straightforward to cali- 
brate and use these virial estimators, one must bear in mind 
the large uncertainties (> 0.4 dex) associated with these es- 
timates and the systematics involved in the calibration and 
usage, which will potentially lead to s ignificant biases of 
these BH mass estimates (e . g.,ICollin et a l. 2006; Shen et al] 
l2008b; Marconi et al. 2008t iDenney et al.ll2009t IKelly et aP 
12009; Shen & Kellv 2010.) . 

The viri al BH mass calib ra tions used in this paper 
are from 'McLure & DunlopI (|2004 H/3 and Mgii), 
Vestergaard & Peterson (200i H/3 and Civ), and 
Vestergaard & O smeil (I2009t Mgii). These calibrations 



have parameters: 








{a,b) = 


(0.672,0.61), 


MD04; H^ 


(3) 


ia,b) = 


(0.505,0.62), 


MD04; Mgii 


(4) 


{a,b) = 


(0.910,0.50), 


VP06; H/3 


(5) 


ia,b) = 


(0.660,0.53), 


VP06; Civ 


(6) 


{a,b) = 


(0.860,0.50), 


VO09; Mgii 


(7) 



In using each of these relations we choose the proper FWHM 
definition adopted in these calibrations. In order to uti- 
lize our new Mgii FWHM measurements (e.g., multiple- 
Gaussian fits with narrow line subtraction, see 33.31 ), we adopt 
the b = 0.62 slope i n the BLR radius-luminosity relation in 
iMcLure & DunlopI (|2004|) . and recalibrate the coefficient a 
such that the Mgii-based estimates are consistent with the H/3- 
based (VP06) estimates on average. We choose this particu- 
lar slo pe b because it was re-calibrated in iMcLure & DunlopI 
using a subsample of reverberation mapping AGNs 
that occupy the high-luminosity regime in the local RM AGN 
sample, arguably better than using the whole RM sample. We 
choose the VP06 H/3 formula to calibrate our Mgii formula 
because the FWHMs of broad H/3 and Mgii were measured 
in a similar fashion (as opposed to the single-Gaussian or 
Lorentzian profile adopted in McLure & Dunlop 2004). This 
new Mgii calibration is 



(«,/?) = (0.740,0.62), S10;Mgii 



(8) 



We note that in a few extremely radio-loud quasars, the continuum lu- 
minosity is significantly boosted by the optical emission from the jet, which 
will then lead to overestimation of the BLR size and the virial BH mass (e.g., 
IWu et al. 2004). The fraction of such objects in our sample is tiny and hence 
we neglect this detail. But we caution the usage of cataloged virial BH masses 
for such objects. 
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In some cases a particular line may be unavailable in the 
spectrum, or the continuum is too faint to measure. Below 
we provide several alternative empirical recipes to estimate a 
virial BH mass. These recipes are calibrated using correla- 
tions among continuum and emission line properties, and are 
only valid in the average sense. We recommend to use them 
only when t he above estimators are unavailable. 

Following l&eene & Hoi (l2005bh . a virial BH mass can be 
estimated based on the FWHM and luminosity of the broad 
Ha line. We found a correlation between the FWHMs (using 
multiple-Gaussian fits) of the broad Ha and H/3 lines simi- 
lar to that found i n previous work (e.g., Gree ne & Holl2005bl 
IShen et aljr2008ah . with the broad H/3 FWHM systematically 
larger than the broad Ha FWHM: 

/FWHMHfl\ 

,(,„5±0.01)]og(5^),9) 

where the slope and intercept are determined using the BCES 
bisector linear regression estimator (e.g., Akritas & Ber- 
shady 1996) for a sample of ^ 2400 quasars with both Ha 
and H/3 FWHM measurements and FWHM errors less than 
500kms~'. Our continuum luminosities at 5100 A have a nar- 
row dynamical range and suffer from host contamination at 
the low-luminosity end, hence instead of fitting a new re- 
lation using our measurements, we adopt the relation be- 
t ween 5 100 A continuum luminosity and Ha line luminosity 
in lGreene & Hol(l2005bt their eqn. 1). The virial mass estima- 
tor based on Ha therefore reads 



log 



Mr 



Mr: 



0.379 + 0.43 log 



Hq 



+2. Hog 



1042 erg s-i 

/FWHMHa 

\ kms"' 



(10) 



where Lhq is the total Ha line luminosity. For quasars in 
our sample, Eqn. (fTOl i yields virial BH masses consistent with 
the VP06 H/3 results, with a mean offset ~ 0.08 dex and a 
dispersion ~ 0.18 dex. 

Similarly, we can substitute the continuum luminosity in 
the above recipes for H/3, Mgii and Civ with the luminosity of 
the particular line used, given that for broad line quasars the 
line luminosity correlates with the continuum luminosity to 
some extent. We determine these correlations using subsam- 
ples of quasars for which both luminosities were measured 
with an uncertainty < 0.03 dex. For the same reason as Ha, 
we do not fit a new relation between the H/3 line luminosity 
an d continuum l uminosity at 5100 A, and we refer to eqn. (2) 
of iGreene & Hq (2005b ) for such a relation. The following 
relations are determined again using the BCES linear regres- 
sion estimator (for both the bisector fit and the Y \X fit, where 
the latter refers to "predict 7 as a function of X"), where the 
line luminosity refers to the total line luminosity'^: 



log ( ) =(2.22 ±0.09) 



ergs 



+ (0.909 ±0.002) log 



ergs"' 



(11) 



bisector ; 



For Mgll, the difference between the total Hne luminosity and the broad 
line luminosity is small enough such that it essentially makes no difference in 
the linear regression results when we use the broad line luminosity instead. 



log( )=(1.22±0.11) 



ergs 



+ (1.016±0.003)log 



-^Mgll 



ergs- 



(12) 
(Y|X) , 



where the scatter of this correlation is ~ 0.15 (0.16) dex for 
^ 44,000 quasars, and 



log ( ) =(4.42 ±0.27) 



ergs 



+ (0.872 ±0.006) log 



log ( -^^^ I =(7.66 ±0.41) 



ergs 



+ (0.863 ±0.009) log 



ergs"' 



Lew 
ergs"' 



(13) 
bisector ; 

(14) 
(Y|X) , 



where the scatter of this correlation is 0.18 (0.2) dex for 
^ 10,000 quasars. Using our measurements, one can also 
estimate these correlations with other linear regression algo- 
rithms. If we substitute logL3()oo and logL^so using Eqns. 
(12) and (14) in the Mgii and Civ estimators, we obtain virial 
BH masses consistent with the original recipes with no mean 
offset (< 0.02 dex) and negligible scatter 0.1 dex). 

There are systematic differences among different versions 
of virial calibrations. For instance, the calibrations for H/3 and 
Mgii in iMcLure & DunlopI (l2004t) u sed the old RM masses 
and vi rial coefficient, while those in IVestergaard & Petersoiil 
('2006') and' Vestergaard & Osme r* (2009')^ used the up dated RM 
masses and virial coefficient (Onken e t al.l |2004|) . Signifi- 
cant uncertainties of the virial coeffi cient still remain (e.g., 
IWoo et all I20T0I: iGraham et alJ 1201 Ih . Moreover, different 
versions of virial calibration for the same line have differ- 
ent dependence on luminosity, and they usually measure the 
line FWHM differently (even though occasionally different 
approaches to measure the FWHM yield the same value dur- 
ing the multi-parameter fits), or prefer an alter native proxy 
(e.g., l ine dispersion) for th e virial velocity (e.g.. lCollin et alj 
I2006t iRafiee &Haiil 12010 ). Currently there is no consen- 
sus on which version of calibration is better It is impor- 
tant to explore various systematics with RM AGN samples 
and statistical quasar samples to determine which is the best 
approach to estimate quasar BH masses with the virial tech- 
nique, and t his is wo rk in prog ress (e.g.. | Onken & Kollmeiei] 
2008; DenneiLeLalJ[2009; Wa ng et al.N2009bl: iRafiee & HaH 
2010). In the mean time, there is strong need to increase the 
sample size and representativity of AGNs with RM measure- 
ments, which anchor these single-epoch virial estimators. 

Here we simply settle on a fiducial virial mass estimate: 
we use H/3 (VP06) estimates for z < 0.7, Mgii (SIO) es- 
timates for 0.7 < z < 1.9 and Civ (VP06) estimates for 
z> 1.9. Fig. [TOl shows the comparison between these virial 
estimates between two lines for the subset of quasars for 
which both line estimates are available and the median line 
S/N> 6. There is negligible mean offset (< 0.01 dex) be- 
tween these virial estimates (half by design), which moti- 
vated o ur choice o f these three calibrations. However, as 
noted in lShen et alj (i2008b), there is a strong trend of decreas- 
ing the ratio of log(M^j"/-'^BH') ^f^h increasing Civ-Mgii 
blues hifts, indicating a po ssible non-virialized component in 
Civ. lRichards et ai](l201 ll) further developed a unified picture 
in which the Civ line has both a non-virial wind component 
and a traditional virial component, and it is plausible that the 
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Fig. 10. — Comparison of virial masses between two different line estimators for the subset of quasars in our sample for which both line estimates are available 
and the median line S/N per pixel > 6. The left panels are one-to-one plots, where the contours are local point density contours estimated with a grid size of 
A = 0.1 on both axes. The right panels show the distribution of mass ratios between two lines, and the mean and Icr from a Gaussian fit to the distribution are 
indicated in the top-left comer. Our new Mgll recipe gives consistent virial mass estimates as the VP06 H/3 recipe (by design), and it also yields consistent virial 
mass estimates as CIV (VP06). These comparison results are similar to those in Shen et al. (2008b, fig. 6), where we compared between H/3 (MD04), Mgll 
(MD04) and CIV (VP06). 



BH mass scaling relation based on Civ is only relevant for 
objects dominated by the virial component. 

3.8. The spectral catalog 

We have tabulated all the measured quantities from the 
spectral fitting in the online catalog of this paper, along with 
other properties'**. The current compilation extends o ur ear- 
lier SP SS Data Release 5 (DR5) compilation (Sh en et all 
l2008bl) by including the post-DR5 quasars, as well as mea- 
surements based on new multiple-Gaussian fits to the lines (as 
discussed above). The format of the catalog is described in Ta- 
bl e 1. Objects are in the s ame order as the DR7 quasar catalog 
in [Schneider et alj ( l2010i) . Below we describe the specifics 



Note that in lShen et alj )2008bl) we only repoiled high-quality measure- 
ments with median S/N> 6 and a reduced < 5 for a single-Gaussian fit to 
the line; here we retain all measurements for completeness. 



of the cataloged quantities. The SDSS terminology can be 
found on the SDSS website'^. Flux measurements were cor- 
rected neither for intrinsic extinction/reddening, nor for host 
contamination. We only report FWHM and velocity shift val- 
ues for detectable line components (i.e., the fitted line flux is 
non-zero), except for cases where the line FWHM and veloc- 
ity offset can be inferred from other lines (such as the narrow 
H/3 line, whose FWHM and velocity offset are tied to those 
of the [Oiii] doublet). 

1. SDSS DR7 designation: hhmmss.ss+ddmmss.s 
(J2000.0; truncated coordinates) 

2-4. RA and DEC (in decimal degrees, J2000.0), redshift. 
Here the redshifts are taken from the DR7 quasar cata- 

http://www.sdss.org/dr7/ 
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log dSchneider etalj|2010h . iHewett & WildI (1201 Ol) pro- 
vided improved redshifts for SDSS quasars. These im- 
proved redshifts are particularly useful for generating 
coadded spectra, but the cataloged DR7 redshifts are 
fine for most of the purposes considered here. 

5-7. Spectroscopic plate, fiber and Modified Julian date 
(MJD): the combination of plate-fiber-MJD locates 
a particular spectroscopic observation in SDSS. The 
same object can be observed more than once with dif- 
ferent plate-fiber-MJD combinations either on a re- 
peated plate (same plate and fiber numbers but differ- 
ent MJD number), or on different plates. The DR7 
quasar catalog typically lists the spectroscopic obser- 
vation with the highest S/N. 

8. TARGET_FLAG_TARGET: the target selection flag 
(TARGET version). 

9. Nspec- number of spectroscopic observations. While we 
only used the default spectrum in our spectral fitting, 
this flag indicates if there are multiple spectroscopic ob- 
servations for each object. 

10. Uniform flag. O=not in the uniform sample'*; l=uni- 
formly selected using t he target selection algorithm in 
[Richards et"an (l2002al) . and flux hmited to i = 19.1 
at z < 2.9 or i = 20.2 at z > 2.9; 2=selected by the 
QSO HIZ branch only in the uniform target selec- 
tion (Richards et al."2002a') and with measured spectro- 
scopic redshift z < 2.9 and / > 19.1. Objects with uni- 
form flag=2 are selected by the uniform quasar target 
algorithm, but should not be included in statistical stud- 
ies; the fraction of such uniform objects is low (< 1%; 
red dots in Fig.[T]l. 

11. M,(z = 2): absolute /-band magnitude in t he current cos- 
mology, /T-corrected to z = 2 following'^ lRichards et alJ 
(2006b). 

12-13. Bolometric luminosity Lboi and its error: computed 
from L5100 (z < 0.7), L3000 (0.7 < z < 1.9), L1350 (z > 
1 .9) using the spectral fits and bolometric corrections'^ 
BC5100 = 9.26, BC3oo q = 5.15 and BCbsq = 3.81 from 
the composite SED in lRichards et alJ (l2006ah . 

14. BAL flag: 0=nonBALQSO or no wavelength cover- 
age; l=Civ HiBALQSO; 2=Mgii LoBALQSO; 3=both 
1 and 2. The LoBALQSO selection is very incomplete 
as discussed in ^ 

For more details regarding the uniform sample selection and its sky 
coverage, see, e.g., Richards et al. 1 2002a, 2006b); Shen et al. (2007). 

" The ^^-corrections h ere include both continuum A'-correction and 
emission-line A'-corre ction jRichards et al. 2006b); while the cataloged ab- 
solute magnitudes in Sch neider et alj (201(31) were A'-corrected for continuum 
only, assuming a power-law continuum. 

The SEDs fo r individual quasars show significant scatter (e.g., 
[Richards et al.l2006afl . so the adopted bolometric corrections are only appro- 
priate in the average sense. Some authors suggest to remove the IR bump in 
the SED in estimating the bolometric corrections (e.g., Marconi et al. 2004), 
where the IR radiation is assumed to come from the reprocessed UV radia- 
tion. This will generally reduce the bolometric corrections by about one third. 
However, we are not correcting for intrinsic extinction of the flux, and we did 
not subtract emission line flux in the Richards et al. composite SED when 
estimating the bolometric corrections, hence using our fiducial bolometric 
connections will not overestimate the bolometric luminosity significantly. 



15. FIRST radio flag: -l=not in FIRST footprint; 0=F1RST 
undetect ed; l=core domina nt; 2=lobe dominant (for de- 
tails, see 'Jiang et al.ll20O7 l, and discussion in Note 
these two classes of radio morphology do not necessar- 
ily correspond to the FR 1 and FR 11 types (see Lin et 
al. 2010 for more details). 

16-17. Observed radio flux density at rest-frame 6 cm f^^m and 
optical flux density at rest- frame 2500 A /2500. Both are 
in the observed frame. 

18. Radio loudness /? = /6cm//2500- 

19-24. L5100, ^-1350 and their errors: continuum luminos- 
ity at 5100 A, 3000 A and 1350 A, measured from the 
spectral fits. No correction for host contamination is 
made (see discussion in 93.61 ). 

25-30. Line luminosity, FWHM, restframe equivalent width 
and their errors for the broad Ha component. 

31-36. Line luminosity, FWHM, restframe equivalent width 
and their errors for the narrow Ha component. 

37-40. Line luminosity, restframe equivalent width and their 
errors for narrow [Nil] A6584. 

41-44. Line luminosity, restframe equivalent width and their 
errors for narrow [S 11] A6717. 

45-48. Line luminosity, restframe equivalent width and their 
errors for narrow [S 11] A673 1 . 

49-50. 6000-6500 A iron restframe equivalent width and its 
error. 

51-52. Power-law slope a\ and its error for the continuum fit 
for Ha. 

53-54. Number of good pixels and median S/N per pixel for 
the Ha region (6400-6765 A). 

55. Reduced for the Ha line fit; -1 if not fitted. 

56-61. Line luminosity, FWHM, restframe equivalent width 
and their errors for the broad H/3 component. 

62-67. Line luminosity, FWHM, restframe equivalent width 
and their errors for the narrow H/3 component. 

68. FWHM of broad H/3 using a single Gaussian fit 
(IShen etal.ll2008H) . 

69-72. Line luminosity, restframe equivalent width and their 
errors for [O III] A4959. 

73-76. Line luminosity, restframe equivalent width and their 
errors for [O III] A5007. 

77-78. 4435-4685 A iron restframe equivalent width and its 
error. 

79-80. Power-law slope a\ and its error for the continuum fit 

for np. 

81-82. Number of good pixels and median S/N per pixel for 
the np region (4750-4950 A). 

83. Reduced for the H/3 line fit; -1 if not fitted. 
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84-89. Line luminosity, FWHM, restframe equivalent width 
and their errors for the whole Mgii profile. 

90-95. Line luminosity, FWHM, restframe equivalent width 
and their errors for the broad Mgii profile. 

96. FWHM of broad Mgii using a single Gaussian fit 
(Shen et al. 2008b). 

97-98. 2200-3090 A iron restframe equivalent width and its 



99-100. Power-law slope ax and its error for the continuum fit 
for Mgii. 

101-102. Number of good pixels and median S/N per pixel for 
the Mgii region (2700-2900 A). 

103. Reduced for the Mgii Hne fit; -1 if not fitted. 

104-109. Line luminosity, FWHM, restframe equivalent width 
and their errors for the whole Civ profile. 

110-111. Power-law slope a\ and its error for the continuum fit 
for Civ. 

112-113. Number of good pixels and median S/N per pixel for 
the Civ region (1500-1600 A). 

1 14. Reduced for the Civ line fit; -1 if not fitted. 

1 15-126. Velocity shifts (and th eir errors) relative to th e systemic 
redshift (cataloged in [Schneider et al.ll2010l) for broad 
Ha, narrow Ha, broad H/3, narrow H/3, broad Mgii, 
and Civ. The velocity shifts for the broad lines are mea- 
sured from the peak of the multiple-Gaussian model fit 
to the broad component'^. Recall that the velocity shifts 
of narrow lines were tied together during spectral fits. 
These velocity shifts can be used to compute the rel- 
ative velocity offsets between two lines for the same 
object, such as the Civ-Mgii blueshift, but should not 
be interpreted as the velocity shifts from the restframe 
of the host galaxy due to uncertainties in the systemic 
redshift. Positive values indicate blueshift and negative 
values indicate redshift; value of 3 x 10^ indicates an 
unmeasurable quantity. 

127-138. Virial BH masses using calibrations of H/3 (MD04), H/3 
(VP06), Mgii (MD04), Mgii (VO09), Mgii (SIO) and 
Civ (VP06). The definitions of the acronym names of 
each calibration can be found in ^3.71 Zero value indi- 
cates an unmeasurable quantity. We use FWHMs from 
a single Gaussian fit to the broad component for H/3 
(MD04) and Mgii (MD04); FWHMs from the multiple- 
Gaussian fit to the broad H^ for H/3 (VP06); FWHMs 
from the multiple-Gaussian fit to the entire Mgii and 
Civ lines for Mgii (VO09) and Civ (VP06) respec- 
tively; FWHMs from the multiple-Gaussian fit to the 
broad Mgii line for Mgii (SIO). See ^for details. 

139. The adopted fiducial virial BH mass if more th an o ne 
estimate is available. See detailed discussion in ^3.71 



SDSS J1 61 742.53+322234.3 
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" The velocity shifts of the broad lines measured from the centroid of a 
single Gaussian fit to the line on average are consistent with those using the 
peak of the multiple-Gaussian fit with negligible mean offset, but they can 
differ (typically by < 200 kms~' ) for individual objects. 



1400 1600 1800 2000 

Rest Wavelength (A) 

Fig. 1 1 . — Two examples of disk emitters with high confidence (bit#0=l). 
Upper: a Balmer disk emitter at z = 0.15. Bottom: a possible CIV disk emitter 
at z= 1.94. 

140. The measurement uncertainty of the adopted fiducial 
virial BH mass, propagated from the measurement un- 
certainties of continuum luminosity and FWHM. Note 
that this uncertainty includes neither the statistical un- 
certainty (> 0.3-0.4 dex) from virial mass calibrations, 
nor the systematic uncertainties with these virial BH 
masses. 

141. Eddington ratio computed using the fiducial virial BH 
mass. 

142. Special interest flag. This is a binary flag: bit#0 
set=disk emitters with high confidence (the vast ma- 
jority are selected based on the Balmer lines); bit#l 
set=disk emitter candidates; bit#2 set=double-peaked 
[O III] AA4959, 5007 lines. These flags were set upon 
visual inspection of all z < 0.89 quasars in the cata- 
log. In particular, disk emitter candidates (bit#l=l) 
are those with asymmetric broad Balmer line profile or 
systematic velocity shifts from the narrow lines; while 
those with high confidence (bit#0=l) show unambigu- 
ous double-peaked (or highly asymmetric) broad line 
profile or large velocity offsets between the broad and 
narrow lines. Fig.fTTIshows two examples of disk emit- 
ters with high confidence. We call these objects "disk 
emitters" even though some of them may be explained 
by alternative scenario s, such as a close S MBH binary 
(e.g., see discussion in lShen & Loebll2010h . 



4. APPLICATIONS 

The spectral measurements described above can be used to 
study the statistical properties of broad line quasars. Here we 
briefly discuss some applications of this spectral catalog. 

4.1. Correlations between emission line properties 

One great virtue of the SDSS DR7 quasar survey is that it 
provides unprecedented statistics for broad-Une quasar prop- 
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Fig. 12. — Statistical properties for H/3 based on our spectral measurements (for all quasars). Contours are local point density contours estimated with a grid 
size of A = 0.1 on both axes (A = 2 on the median S/N axis). Contour levels are equally spaced linearly with their values marked for the outermost several 
contours. The strong an ti-cor relation between the power-law continuum slope and luminosity reflects the increasing host galaxy contamination towards fainter 
quasar luminosites (see j|3.6l and Fig.|9). 
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Fig. 13. — Statistical properties for Mgll based on our spectral measurements (for all quasars). Contours are local point density contours estimated with a 
grid size of A = 0.1 on both axes (A = 2 on the median S/N axis). Contour levels are e qually spaced l i nearly with their va lues marked for the outermost several 
contours. The re are several know n correlations, such as the Mgll Baldwin effect (e.g.. lBaldwinlll97"7tlCroom et alj|2002l) . and the correlation between FWHM 
and EW (e.g.. lDong et aU2009al) . 
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Fig. 14. — Statistical properties for CIV based on our spectral measurements (for all quasars). Contours are local point density contours estimated with a 
grid size of A = 0.1 on both axes (A = 2 on the median S/N axis). Contour levels are equally spaced hnearly with their values marked for the out ermost several 
contours. Ther e are several coiTelations involving EW, FWHM and luminosity that may be different manifestations of the same phenomenon (e.g.. BaldwinI 1 9771 : 
[Richards et al.l 2002b). The clustering of a population of low S/N objects is caused by the quasar target selection, i.e., quasars are targeted to a fainter limiting 
magnitude at z > 2.9 (see Fig.[T). 
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Fig . 15. — Distrib ution of quasars in the projected 4DE1 parameter space 
(e.g.. ISulentic et aT] i20Q0). Contours are local point density contours esti- 
mated with a grid size of A.v/Ay = 0.1/500 (black) and A, /A,. = 0.2/800 
(red), where the black contours are for all quasars and the red contours are 
for radio-loud (R > 10) quasars only. 

erties. To demonstrate this, Figs. [T2lfT4l show some statisti- 
cal properties of quasars using our spectral measurements for 
H/3, Mgii, and Civ, respectively. These figures show the typ- 
ical values of these properties for SDSS quasars as a quick 
reference. We do not show a similar figure for Ha because 
quasars with Ha coverage represent only a tiny fraction of the 
whole sample, and because host contamination is more severe 
for these low-redshift quasars. 

There are correlations among the properties shown in Figs. 
n"2lfT4l Some of these correlations are not due to sel ection ef- 
fects. For instance, the well-known Baldwin effect (iBaldwinI 
Il977l) . i.e., the anti-correlation between line EW and contin- 
uum luminosity, is clearly seen for Civ and Mgii. There are 
also strong correlations between EW and FWHM for Mgii 
(e.g., Dong et al. 2009b) and Civ, which are not due to any 
apparent selection effects. The statistics of our catalog now 
allows in-depth investigations of these correlations when bin- 
ning in different quantities such as redshift or luminosity, and 
to probe the origins of these correlations. However, there are 
some apparent correlations which are likely due to selection 
effects inherent in a flux-limited sample, or host contami- 
nation. For instance, the apparent anti-correlation between 
ax and logLsioo, and the mild negative Baldwin effect below 
logLsioo < 45 for H/3 seen in Fig. [12] are most likely caused 
by increasing host contamination towards fainter luminosities 
(see ^3.61 ). Moreover, the spectral quality (mainly S/N) has 
important effects on the measured quantities and may bias the 
measurements at the low S/N end. Thus one must take these 
issues into account when using the catalog to study correla- 
tions among various properties. The detailed investigations of 
various correlations will be presented elsewhere. 

Fig. [15] shows the so-called 4DE1 projection in the H/3 
FWHM versus R^i^n = E \ yFBn44^-4fiS4 /EWh /9 space (e.g., 
ISulentic et al] l200d 120021 ; IZamfir et al.l l20ld ). which is an 



extension of the eigenvector space for quasar properties sug- 
gested by ' Boroson & Greenl(ll992h . Objects with H/3 FWHM 
> 4000 kms"' (i.e., population "B" in the terminology of the 
4DE1 parameter space of Sulentic and collaborators) have a 
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Fig. 16. — Correlation between L[oill]5007 and Z.5100. Contours are lo- 
cal point density contours estimated with a grid size of A = 0.1 on both 
axes. The dashed line is the linear regression fit treating logLsioo as the 
independent variable; the best-fit relation is log(L[oin]5007/erg s"' ) = 7.76 + 
0.771og(L5ioo/ergs"') and the scatter around this relation is ~ 0.38 dex. 
The dash-dotted line is the bisector linear regression fit; the best-fit rela- 
tion is log(L[oin]5007/ergs"') = -17.54 + 1.341og(L5ioo/ergs"') and the scat- 
ter around this rel ation is ~ 0.4 dex. The solid line is the mean Hnear relation 
described in Eqn. )15t . 

tendency to have weaker relative iron emission strength for 
larger FWHMs. The black contours show the distribution 
of all quasars while the red contours show the distribution 
of radio-loud (R > 10) quasars. It appears that the radio- 
loud contours are more vertically elongated, broadly consis- 
tent with the phenomenological c lassification scheme based 
on the 4DE1 param eter space (e.g.. ISulentic et alJ2000l 120021 : 
IZamfir et alJl2010D . The physics driving these characteristics 
in the parameter space is currently not clear, and deserves fur- 
ther study. 

Fig.[T6lshows the correlation between the [Oiii] A5007 lu- 
minosity and the continuum luminosity at 5100 A. This cor- 
relation is usually used to estimate the bolometric luminos- 
ity using the [O iii] A5007 luminosity as a sur rogate for type 
2 quasars (e.g, Kauffmann et al. 2003; Zaka mska et alj|2003l ; 
iHeckman etd]r2004: .Reyes et al.,.2008.) . While the correla- 
tion is a pparent, it has a large scatter, as noted in earlier stud- 
ies (e.g.. IHeckman et al.ll2004t iReves et al.ll2008l) . The mean 
linear relation is: 



logZ'PIIIlASOO? ~ logL51()0-2.5 



(15) 



with a s catter ~ 0.35 dex. U sing the bolometric correc- 
tion from lRichards et al] (|2006a), a crude conversion between 
i'[Oin]A5007 and the quasar bolometric luminosity is; Lboi ~ 
3200L[om]A5007- 

4.2. Emission line shifts 

Fig. [12] shows the distributions of velocity shifts between 
various emission lines. Recall that the velocity of the broad 
lines is measured from the peak of the multiple-Gaussian 
fit. The left panel of Fig. [17] shows the distributions of ve- 
locity shifts between the broad Balmer lines and the narrow 
lines. The means of these distributions are consistent with 
zero, hence there is no offset in the mean between the broad 
and narrow Balmer lines (cf.. iBonning et al]|2007h . We note 
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Fig. 17. — Velocity shifts between different pairs of lines. Positive values indicate blueshifts. 



that if we did not account for the blue wings of the nar- 
row [O III] AA4959,5007 lines during spectral fitting, there 
would be a net redshift of the order of ^ 100 kms~' be- 
tween the broad H/? line and [Oiii], which is inconsistent 
with the results for Ha versus [Sii]. The right panel of Fig. 
[TtI shows the velocity offsets between Mgii and [O iii] and 
between Civ and Mgii. The Mgii line shows no mean offset 
from [Oiii], while the Civ line shows a systematic blueshift 
of ^ 600 kms"' wit h respect to Mgii ( e.g., Gaskell 1982; 
l^tler & Fan 1992; Richards et al.ll2002bh : see lRichards etal.1 
( 1201 1) for further discussion. 

It is interesting to note that many of the objects in the 
wings of the velocity offset distributions of the broad Balmer 
lines vs the narrow lines are either strong disk - emitters (e.g., 
Chen et all [T989t lEracleous & HalpernI 119941; IStrateva et all 
2003h . or have the broad component systematically offset 
from the narrow line center; in other cases the apparent large 
shifts were caused by poor fits to noisy spectra. 

5. SUMMARY 

We have constructed an SDSS DRV quasar catalog in which 
we tabulate various properties. In this catalog we compiled 
continuum and emission line properties for Ha, H/3, Mgii, 
and Civ based on our spectral fits. We also included radio 
properties, and flagged quasars of special interest, such as 
broad absorption line quasars and disk emitters. We also com- 
piled virial BH mass estimates using these spectral measure- 
ments. This catalog can be used to study correlations among 
properties of optically selected quasars, and the active black 
hole mass function in quasars (Shen et al., in preparation). 
We performed various tests and found that our automatic fit- 
ting procedure to emission lines performed reasonably well. 
However, one must take into account the possible effects of 
selection and S/N, as well as the systematics involved in con- 
verting the measured quantities to derived quantities, when 
using these measurements to study quasar properties. In par- 
ticular, we do not encourage direct interpretations based on 
derived quantities (such as virial BH masses and bolometric 
luminosity) without accounting for the difference between the 
estimated value and the true value for these quantities. Such 



direct interpretations will usually lead to biased or even spu- 
rious results. 

We make this catalog publicly available online^°, where 
we also provide supplemental materials (such as dereddened 
spectra, QA plots, etc) and important future updates of this 
compilation. 
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PROPERTIES OF SDSS DR7 QUASARS 

TABLE 1 FITS Catalog Format 



Column Format Description 

1 STRING SDSS DR7 designation hhmmss.ss+ddmmss.s (J200().0) 

2 DOUBLE Right ascension in decimal degrees (J2000.0) 

3 DOUBLE Declination in decimal degrees (J2000.0) 

4 DOUBLE Redshift 

5 LONG Spectroscopic plate number 

6 LONG Spectroscopic fiber number 

7 LONG MJD of spectroscopic observation 

8 LONG Target selection flag (TARGET version) 

9 LONG Number of spectroscopic observations 

10 LONG Uniform selection flag 

11 DOUBLE M,(z = 2) [ft = 0.7. f^o = 0.3, J^a = 0.7, /^-corrected to z = 2, following Richards et al. (2006b)] 

12 DOUBLE Bolometric luminosity [log(Lboi /ergs"')] 

13 DOUBLE Uncertainty in log Lboi 

14 LONG BAL flag {0=nonBALQSO or no wavelength coverage; 1=CIV BALQSO; 2=MgII BALQSO; 3=bofli 1 and 2) 

15 LONG FIRST radio flag (-l=not in FIRST footprint; 0=FIRST undetected; l=core-dominant; 2=lobe-dominant) 

16 DOUBLE Observed radio flux density at rest-frame 6 cm /y.6cm [mJy] 

17 DOUBLE Observed optical flux density at rest-frame 2500A [log(/,^ 250o/ergs"'cm"^Hz"')] 

18 DOUBLE Radio loudness R = /^ gc,^//^ 2500 

19 DOUBLE Monochromatic luminosity at 5100A [log(L5ioo/ergs"')] 

20 DOUBLE Uncertainty in log L5 

21 DOUBLE Monochromatic luminosity at 3OOOA [log(L3ooo /ergs"')] 

22 DOUBLE Uncertainty in logLsooo 

23 DOUBLE Monochi'omatic luminosity at 1350A [log(Li35o/ergs"' )] 

24 DOUBLE Uncertainty in logLisso 

25 DOUBLE Line luminosity of broad Ha [log(L/ergs"')] 

26 DOUBLE Uncertainty in log Lhq. broad 

27 DOUBLE FWHM of broad Ha (kms"' ) 

28 DOUBLE Uncertainty in the broad Ha FWHM 

29 DOUBLE Restframe equivalent width of broad Ha (A) 

30 DOUBLE Uncertainty in EWhq .broad 

31 DOUBLE Line luminosity of nan'ow Ha [log(L/ergs"')] 

32 DOUBLE Uncertainty in log Lhq, narrow 

33 DOUBLE FWHM of narrow Ho (kms"' ) 

34 DOUBLE Uncertainty in the naiTow Ha FWHM 

35 DOUBLE Restframe equivalent width of narrow Ha (A) 

36 DOUBLE Uncertainty in EWna. narrow 

37 DOUBLE Line luminosity of [N II] A6584 [logCL/ergs"' )] 

38 DOUBLE Uncertainty in logL[f,jii]5584 

39 DOUBLE Restframe equivalent width of [N II] A6584 (A) 

40 DOUBLE Uncertainty in EWiniij6584 

41 DOUBLE Line luminosity of [S II] A6717 [log(L/ergs-' )] 

42 DOUBLE Uncertainty in logL[sii]S7i7 

43 DOUBLE Restframe equivalent width of [S II] A6717 (A) 

44 DOUBLE Uncertainty in EW[sii]S7 17 

45 DOUBLE Line luminosity of [S II] A6731 [log(L/ergs-' )] 

46 DOUBLE Uncertainty in logL[sii]673i 

47 DOUBLE Restframe equivalent width of [S II] A673 1 (A) 

48 DOUBLE Uncertainty in EW[siij673i 

49 DOUBLE Restframe equivalent width of Fe within 6000-6500A (A) 

50 DOUBLE Uncertainty in EWfc.Hq 

51 DOUBLE Power-law slope for the continuum fit for Ha 

52 DOUBLE Uncertainty in ajja 

53 LONG Number of good pixels for the restframe 6400-6765A region 

54 DOUBLE Median S/N per pixel for the restframe 6400-6765A region 

55 DOUBLE Reduced for the Ha line fit; - 1 if not fitted 

56 DOUBLE Line luminosity of broad H^ [log(L/ergs-' )] 

57 DOUBLE Uncertainty in log Lh/9 .broad 

58 DOUBLE FWHM of broad H/3 (kms"') 

59 DOUBLE Uncertainty in the broad H/3 FWHM 

60 DOUBLE Restframe equivalent width of broad H/3 (A) 

61 DOUBLE Uncertainty in EWh/3, broad 

62 DOUBLE Line luminosity of narrow H/3 [log(L/ergs"')] 

63 DOUBLE Uncertainty in logLn/j, narrow 

64 DOUBLE FWHM of narrow H/3 (kms"' ) 

65 DOUBLE Uncertainty in the narrow H/3 FWHM 

66 DOUBLE Restframe equivalent width of narrow H/3 (A) 

67 DOUBLE Uncertainty in EWh/3 narrow 

68 DOUBLE FWHM of broad H/3 using a single Gaussian fit (kms"' ) 

69 DOUBLE Line luminosity of [O III] A4959 [log(L/ergs"' )] 

70 DOUBLE Uncertainty in logL[om]4959 

71 DOUBLE Restframe equivalent width of [O III] A4959 (A) 

72 DOUBLE Uncertainty in EW[oin]4959 

Continued on Next Page. . . 
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SHEN ET AL. 
TABLE 1 - Continued 



Column 


Format 


Description 


73 


DOUBLE 


Line luminosity of [0 III] A5007 [log(L/ergs~')] 


74 


DOUBLE 


Uncertainty in logZ,[om]5007 


75 


DOUBLE 


Restframe equivalent width of [O III] A5007 (A) 


76 


DOUBLE 


Uncertainty in EW[om]5007 


77 


DOUBLE 


Restframe equivalent width of Fe within 4435-4685A (A) 


78 


DOUBLE 


Uncertainty in EWfcH/S 


79 


DOUBLE 


Power-law slope for the continuum fit for H/9 


80 


DOUBLE 


Uncertainty in oh^ 


81 


LONG 


Number of good pixels for the restframe 4750-4950A region 


82 


DOUBLE 


Median S/N per pixel for the restframe 4750-4950A region 


83 


DOUBLE 


Reduced for the H/3 hne fit; -1 if not fitted 


84 


DOUBLE 


Line luminosity of the whole Mgll [log(L/ergs~^ )] 


85 


DOUBLE 


Uncertainty in logLMgii,whoie 


86 


DOUBLE 


FWHM of the whole Mgll (kms"') 


87 


DOUBLE 


Uncertainty in the whole Mgll FWHM 


88 


DOUBLE 


Restframe equivalent width of the whole Mgll (A) 


89 


DOUBLE 


Uncertainty in EWMgii.whoie 


90 


DOUBLE 


Line luminosity of broad Mgll [log(L/ergs~')] 


91 


DOUBLE 


Uncertainty in IosLmoti hman 


92 


DOUBLE 


FWHM of broad Mgll (kms"') 


93 


DOUBLE 


Uncertainty in the broad Mgll FWHM 


94 


DOUBLE 


Restframe equivalent width of broad Mgll (A) 


95 


DOUBLE 


Uncertainty in EWMgii.broad 


96 


DOUBLE 


IWHM of broad Mgll using a single Gaussian fit (kms~') 


97 


DOUBLE 


Restframe equivalent width of Fe within 2200-3090A (A) 


98 


DOUBLE 


Uncertainty in EWpcMgii 


99 


DOUBLE 


Power-law slope for the continuum fit for Mgll 


100.... 


DOUBLE 


Uncertainty in ccMgii 


101.... 


LONG 


Number of good pixels for the restframe 2700-2900A region 


102 


DOUBLE 


Median S/N per pixel for the restframe 2700-2900A region 


103 


DOUBLE 


Reduced for the Mgll line fit; -1 if not fitted 


104 ... . 


DOUBLE 


Line luminosity of the whole CIV [log(L/ergs^')] 


105 ... . 


DOUBLE 


Uncertainty in logLciv 


106 ... . 


DOUBLE 


FWHM of the whole CIV (kms"' ) 


107.... 


DOUBLE 


Uncertainty in the CIV FWHM 


108.... 


DOUBLE 


Restframe equivalent width of the whole CIV (A) 


109.... 


DOUBLE 


Uncertainty in EWcrv 


110.... 


DOUBLE 


Power-law slope for the continuum fit for CIV 


111.... 


DOUBLE 


Uncertainty in acrv 


112.... 


LONG 


Number of good pixels for the restframe 1500-1600A region 


113 ... . 


DOUBLE 


Median S/N per pixel for the restframe 1500-1600A region 


114.... 


DOUBLE 


Reduced tor the CIV fit; -1 if not fitted 


115 


DOUBLE 


Velocity shift of broad Ho^ (kms ' ); 3d5 if not measurable 


116.... 


DOUBLE 


Uncertainty in Vhq,, broad 


117 


DOUBLE 


Velocity shift of narrow Ho^ (kms^^); 3d5 if not measurable 


118 


DOUBLE 


Uncertainty in Vn^^ n^n-dw 


119 ... . 




vciuciiy siiiii ui Lfiutiti 11/-' juj 11 iiui iiicaoiuaui^ 


120 




uiiLciiaimy iii vp^^f broad 


1 Z 1 .... 




Velocity shift of narrow H/3 (kms '); 3d5 if not measurable 


iZZ .... 




Uncertainty in Vh/3. narrow 


123 


DOUBLE 


Velocity sliirt of broad Mgll (kms ' ); 3a5 it not measurable 


124 .... 


DOUBLE 


Uncertainty in VMgii.broad 


125.... 


DOUBLE 


Velocity shift of CIV (kms"' ); 3d5 if not measurable 


126 


DOUBLE 


Uncertainty in Vciv 


127 ... . 


DOUBLE 


Virial BH mass based on H/3 [MD04, log(MBH.vii /M0)] 


128 


DOUBLE 


Measurement uncertainty in logMaH.vir (H/3> MD04) 


129 


DOUBLE 


Virial BH mass based on H^ [VP06, log(MBH,vii/M0)] 


130 


DOUBLE 


Measurement imcertainty in logMBH.vir (H/3, VP06) 


131.... 


DOUBLE 


Virial BH mass based on Mgll [MD04, iog(MBH,vir/A^©)] 


132 


DOUBLE 


Measurement uncertainty in logMsH.vir (Mgll, MD04) 


133.... 


DOUBLE 


Virial BH mass based on Mgll [VO09^ log(/WBH.vir/M0)] 


134.... 


DOUBLE 


Measurement uncertainty in logMgn vir (Mgll, VO09) 


135.... 


DOUBLE 


Virial BH mass based on Mgll [SIO, log(MBH.vii/M0)] 


136.... 


DOUBLE 


Measurement uncertainty in logMaH vir (Mgll, SIO) 


137.... 


DOUBLE 


Virial BH mass based on CIV [VP06, log(MBH,vir/M0)] 


138.... 


DOUBLE 


Measurement uncertainty in logMBH.vir (CIV, VP06) 


139.... 


DOUBLE 


The adopted fiducial virial BH mass [log(MBH.vir/'W0)] 


140 ... . 


DOUBLE 


Uncertainty in the fiducial virial BH mass (measurement uncertainty only) 


141 ... . 


DOUBLE 


Eddington ratio based on the fiducial virial BH mass [logtLboi /itdd)! 


142. .. . 


LONG 


Special interest flag 


Continued on Next Page. . . 
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Column Format Description 



Note. — (1) Objects are in the same order as in the DR7 quasar catalog ^Sc hneider et alj | 20 1()l); (2) A"— corrections are the same as in [Richards et alJ *2006b|); (3) Bolometric 
luminosities computed using bolometric coiTections in Richards etal. i 2006a) using one of the 5100A, 3000A, or 1350A monochromatic luminosities depending on redshift; (4) 
Uncertainties are measurement errors only; (5) Unless otherwise stated, null value (indicating unmeasurable) is zero for a quantity and —I for its associated enor. 



